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Abstract 

A variational principle for two-fluid mixtures is proposed. The Lagrangian 
is constructed as the difference between the kinetic energy of the mixture and 
a thermodynamic potential conjugated to the internal energy with respect to 
the relative velocity of phases. The equations of motion and a set of Rankine- 
Hugoniot conditions are obtained. It is proved also that the convexity of the 
internal energy guarantees the hyperbolicity of the one-dimensional equations 
of motion linearized at rest. 



1 Variational approach to the description of homogenous two- velocity 
media 

The variational approach to the construction of two-fluid models was used 
by many authors (A.Bedford & D.S. Drumheller (1978), V.L.Bcrdichcvsky 
(1983), J.A. Geurst (1985,1986), H. Gouin (1990)). Here we give its gener- 
alisation for the case of homogeneous two-fluid mixtures. A physical example 
of such flow is a motion of a mixture of two gases with quite different molecular 
weights. The Hamilton's principle is applied for perfect fluid motions. In order 
to obtain the equations of motion of the mixture from this variational princi- 
ple, we neglect the dissipative effects. We will then consider only mechanical 
processes by suppressing thermal evolution. We suppose that the homogeneous 
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mixture motion is well represented by the velocities of its components Ui, Ui , 
the average densites pi, p 2 and the total internal energy U . The total energy 
of a two-velocity medium is written in the form: 

In order to define the internal energy of one- velocity media, a moving coordinate 
system, where the elementary volume of the continuum is at rest, is considered. 
The total energy of the continuum with respect to this system is called the 
internal energy of the motion. For a two- velocity medium, there is no coordinate 
system, within the framework of which any motion could be disregarded. This is 
the reason why the standard definition of internal energy leads to its dependence 
on the relative motion of components. 

Let w = u 2 — vi\ , w = | w | . We propose the following extended form of 
Hamilton's principle of least action for the two- velocity systems: 

SI = , I = J J Jjpi ^f- + P 2 ^f- - W ( Pl , P2 , w))dx dt (1.1) 

with the kinematic contraints 

^+div( Pi Ui)=0. (1.2) 

Here [ ti, t2\ is a time interval , T> is a domain in the physical space, the potential 

W{p\ , p2 ,w) is connected with the internal energy U by the partial Legendre 

transformation with respect to the variable w: 

dW dW 
U(p 1} p 2 ,i) = W(p 1} p 2 ,w) -w-q— = W + i w, with t = ~~g^'- ( 1 ' 3 ) 

The internal energy of the medium U (pi , pi , i) is generally a convex function 
of its variables. This assumption corresponds to the condition of thermodynamic 
stability. As a consequence, W(pi ,p 2 ,w) is a convex function with respect to 
variables p\ , p 2 and a concave function with respect to variable w. A simple 
case is associated with 

W(p! , p 2 , w) = e{pi , p 2 ) , 
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U (pi , P2 , i) = e (pi , p 2 ) + —, r, (1.4) 

l a{pi p2) 

where e{p\ , P2) is a convex function with respect to p\ , p 2 and a {p\ , p 2 ) 
is a positive function. If the relative velocity w is small enough, the energy U 
is convex. The example (1.4) is a reasonable approximation of a general case as 
W is an analytic function of the velocities ttj (and hence, an analytic function 
of w 2 only). The term in equation (1.4) which is quadratic with respect to 
the relative velocity w, can be considered as the energy due to an added mass 
effect. We note that this quadratic dependence is usually used in the theory of 
bubbly liquids (V.L. Berdichevsky, 1983, J. A. Geurst, 1985, 1986). To derive the 
governing equations and the Rankine-Hugoniot conditions, it is not necessary 
to focus on the particular case (1.4). 



2 Governing equations 



We introduce Lagrange coordinates Xi for each component: 
It follows from (2.1) that 

* - - <§>" < f >• (2 2 » 

where the operation A < / > denotes the product of the tensor A by the vector 
/• 

Introducing the Lagrange multipliers tpi(t,x) , ip 2 (t,x) , corresponding to the 
balance of masses (1.2) , we consider the Lagrangian L of the system 

L =M^H 2 - ^r) + /^l^l 2 - ^)-W(p u p 2 ,w)- (2-3) 

dX- 

Thc formulae (2. 2), (2. 3) give the Lagrangian L as a function of variables —^j- , 

dXi d <pi d ifi . ...... 

-—r , — — , —^r, pi,p 2 - Calculating the corresponding variational derivatives, 
ox at ox 

we find the governing equations of motion : 

^ + ^%*-*>+V(^)=0, (2.4) 
dt Ox dpi 
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= - (-1) 4 — -5 , z = l,2 

Pi aw w 

and, of course, the mass conservation laws (1.2). Here and later "*" denotes the 
transposition. For the case (1.4) equations (2.4) have the form : 

dt p l dx \ dpi 2 dpi J 

K i = u i + (-l) w. 

The system (1.2) , (2.4) yields the momentum conservation law and the energy 
conservation law, corresponding to the homogeneity of the Lagrangian with 
respect to space and time variables: 

d ( _ ^ \ ,. , _ ^ _ dW w®w* 

TT7 PlWl + /0 2 M2 + dlV (piUl ® U x + p 2 U 2 ®U 2 ~ — h 

ot\ J aw w 

dW dW \ 
+ (Pi 7 ^+P2 7 —-W)l)=0, (2.5) 
dpi dp 2 J 

£ ^PlH '+2*11*1 ^J+dlV L_L + ^+^2 (4 L + ^J- 

— ^— («2 <8> «2 - ® Ml) < - > =0, (2.6) 
aw w J 

where ® denotes the tensor product and / is the unit tensor. 



The question of the hyperbolicity of the system (1.2), (2. 4) is of great interest. 
Unfortunately, the multi-dimensional case is not very simple and this is why 
we restrict our attention on the one-dimensional case. Moreover, to simplify 
the calculations, we consider the case (1.4) and linearize our system in the 
neighbourhood of the equilibrium state u\ = 0% — , Pi,p 2 - Straightforward 
calculations give the following result: if e (p\ , p 2 ) is a convex function and 
a (pi , p 2 ) > then linearized at rest, the system (1.2), (2.4) is hyperbolic. 

3 The Rankine-Hugoniot conditions 

The non-linearity and hyperbolicity of the system (1.2), (2. 4) imply the ne- 
cessity to obtain Rankine-Hugoniot conditions across shocks. In multi-dimensional 
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case the conservation laws (1.2), (2. 5), (2. 6) are not sufficient to obtain the whole 
set of such relations. In what follows, we show that the additional jump condi- 
tions can be derived from the variational principle. 



Let us define variations of particules deduced from the relation x = 
$i (Xi , t , Si) and its inverse Xi = $j (x , t , £j). Here , i = 1,2 
are small parameters defined in a neighbourhood of zero. One defines virtual 
displacements 5iX and SiXi by 

,5.x = ^-1(^,4,0) , <5 4 X, = ^-1(2, t,0). (3.1) 

— * dx 

It is clear that = —F^ < Six >, where Fj = — — . For any variable 

^ ax, 

aj (a? , t) we define their eulerian perturbations &i (x , t , £j) and lagrangian 
perturbations cifi (^ , i, £j) = ai($>i(Xi, t, £$), i, £j). The variations £jO!i and 
Jittj of Qj are defined by 

koti = -7^- (x,t,0), 6ica = -r^- (Xi,t,0). 

It follows from the above definitions that 

bion = Sidi - < Si x > . (3.2) 

In particular, 

SiPi = -div(p. Six) , Siili — ^ Six - < Six > . (3.3) 

We define 

£i^0 £j 

Then 

Sil= / (^»(^-^— )d2 dt. (3.4) 

It follows from (3.3) (3.4) that 

* ' = £ X <T(| w-§ < * * x) ** 
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/*! JT> 



d - dW \u\ 2 

+-g^(p i S l x* Ki) + div(p l u l (<5 4 x* Ki) + Pl 5ix{— - L ^~)) ) dx dt. 



As a consequence, we obtain equations of motion (2.4) and the jump conditions: 

, (3.5) 



dW \uA 2 

Pi{n* Six) - + Pt (n* Ui) (K* 5 t x) - D n Pi {K*5ix) 



where n is the unit normal vector to the shock surface and D n is the normal 
velocity of the shock. Since 8$ is not continuous across the shock, it is 
not straightforward to obtain Rankinc- Hugoniot conditions in terms of desired 
quantities. Nevertheless, taking into account that across the shock 

[ Pi F* < n >} = 0, 

such difficulties can be overcame. We finally obtain from (3.5): 

, (3.6) 



[ Ki - (K* n)n } = 0. (3.7) 

In the limiting case when the velocities of the components coincide, these condi- 
tions reduce to the conservation of the tangential component of the velocity and 
the conservation of the Bernoulli constant across the shock. We note also that 
conditions (3.6), (3.7) are obtained from the variational principle (1.1) without 
any assumption on the flow properties. But it can be shown directly that they 
correspond to the jump conditions for the additional conservation laws admitted 
by the system (1.2), (2.4) : 

roti? l = 0, ^ + V*(^-l\ Ul \ 2 + KZ Ui) = 0. (3.8) 

The conservation laws (2.5), (2.6) also imply momentum balance and energy 
balance at the shock. Finally note that all possible jump conditions are of great 
interest but the correct choice of the jump conditions depends on the physics of 
the problem. 
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